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Abstract 

In this paper, we derive explicit formulas for the surface averaged first exit time of a discrete 
random walk on a finite lattice. We consider a wide class of random walks and lattices, including 
random walks in a non-trivial potential landscape. We also compute quantities of interest for 
modelling surface reactions and other dynamic processes, such as the residence time in a subvolume, 
the joint residence time of several particles and the number of hits on a reflecting surface. 
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I. INTRODUCTION 



The theory of random walks on lattices is not only a beautiful mathematical object, it 
is a ls0 Us e tul in _ d om ain s of ph y sicsQ, including po t en tW ^ g statistic, 
field theory or biophysics j^J. Another natural application is the diffusion of adatoms 
and vacancies on a crystal surface 0,0, 7|. 

Among the numerous issues involved in the study of such lattice random walks, one im- 
portant area is concerned with random walks on finite lattices. There are two important 
reasons for that special interest. First, true physical systems are not infinite, so that explicit 
boundary conditions have often to be taken into account in order to properly describe situ- 
ations in which confinement can be relevant. Second, exact solvable random walk problems 
in bounded domains are v ery rare, making this theoretical field an important problem in its 
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Recently, Blanco and Fournier [l6| reported an important general result concerning the 
mean first exit time of Pearson random walks jl^ - that is continuous space and time random 
walks, with a given frequency of reorientation A of the direction of the constant velocity v - 
in a bounded domain. They showed that the mean first exit time of a random walk starting 
from the boundary of a finite domain is independant of the frequency A of redirection, and 
is simply related to the ratio of the domain's volume V over the surface S of the domain's 
boundary. The corresponding equation is (in three dimensions): 

(1) 

where v is the speed of the walker. This result was extended by Mazzolo to the higher- 
order moments of the first exit time, and by Benichou et al. Q| to general diffusion processes 
in a non-uniform energy landscape. 

In this paper, we show how these results can be extended to the important case of discrete 
space and time random walks on a finite lattice. In particular, we obtain very simple explicit 
expressions for mean exit times and mean residence times averaged over the surface of the 
considered domain, for rather general random walks. More precisely, the article is organized 
as follows. In section m we define the model under study and the basic averages involved in 
the sequel. The section HTT1 is devoted to the study of first exit time moments. In section HV| 
we generalize this approach and explicitly calculate the mean residence time in a subvolume 
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of the domain. The section IVl presents an analysis of the mean number of hits of a reflecting 
surface. In section IVI1 we generalize the previous results to the very important case of 
a random walk in presence of a potential. Finally, in section IV11[ we consider the joint 
residence time of several particles in subdomains of the lattice and derive results that can 
be applied, for instance to the theory of heterogenous catalysis. 



II. THE MODEL 



Let us start with the definition of our model. First, we have a lattice, which may be of 
any dimension or connectivity: for instance, we can as well apply our results to the cubic 
3-D lattice than to the triangular 2 — D lattice. We study the motion of a random walker 
(with memory: the random walker has a probability of switching direction A each time it 
visits a site; note that we can go back to the model without memory by taking A = 1.) The 
random walker starts from the boundary of a domain (see fig E]) 



FIG. 1: A typical volume: The surface is the dashed line; to the right, the conventions for the 
surface 



The position of the random walker is denoted by r(£); The speed of the random walker is 
v{t). The only values it can take are the difference between the positions of two neighbours. 
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The precise rules of the model thus write: 

f{t + 1) 



= r(t) + v{t) 



(2) 



v(t + 1) 



v(t) with probability 1 — A 
random with probability A 



(3) 



We will also define conditional probabilities: p(r f , v', t\f, v) is the probability for the 
random walker to start from the position r, with the speed v and arrive after a time t at 
the position r* with the speed if. We also note p(f, v, t) the probability for reaching the 
boundary after a time t, when starting from the point r with the speed u; We have 



The sum here means a sum over all the lines which cross the boundary, and to which we 
can associate a point T, in (inside the boundary), a point S out (outside the boundary), and 
a speed v^, which points from the outside to the inside, (see fig^) This convention means 
that we consider that the random walker has reached the boundary as soon as it reaches the 
point inside the boundary, with a speed pointing out. This also means that the time needed 
to reach the boundary is equivalent to the number of (not necessarily distinct) sites visited, 
not including the site the random walker begins on. Thus, with these probabilities, for any 
quantity <p(t) depending on the first exit time t we can define its average on t: 




(4) 



oo 




(5) 



t=o 



For example, we can define the average exit time with simply (p = t 




(6) 



t=0 



Another such average which will be useful is the Laplace transform: 



oo 




(7) 



i=0 



We will then define two spatial averages of tp: the first one is the surface average: 




(8) 
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In particular, < t >s is the mean time needed to return to the surface, or, alternatively, 
the mean number of (not necessarily different) sites visited between entrance and exit. Note 
that the first site, which is out of the volume, is not counted, due to the definition of the 
probabilities. S is the surface, or the number of lines crossing the boundary. 

There is a second useful average which may be defined, the volume average: 

^v = y^J2^v) (9) 

r£V,v 

V is the volume, i.e. the number of sites inside the boundary, and <jd is the coordination 
number of the lattice. For a hypercubic lattice in dimension D, it is simply 2D. 

An example of such a volume average is (t) v , which is the mean number of sites (still 
not necessarily distinct) visited by the random walk (starting from a random point in the 
volume) before exiting, not including the site it starts from. To finish this introduction, let 
us note a few important relations: 

Tp(T, in , -fa) = <p(£ in , -vz, 0) (10) 

This because the random walker starting at this position automatically leaves the volume 
at time 0. We will often have to compute terms like YltLo <p(t)p(r — v, v, t + 1). We can 
notice that, if f £ V, then p(f — v,v,0) = 0. Indeed, p(f,v) is equal to one if r — and 
v = —v-e, and zero else. And, from the definition of S iri and v%, we can see that there is no 
r £ V such that f — v = S in , and v = —v-e. Thus, we have: 

oo oo 

- v,v,t + l) = ^2 <p(t - l)p(r- v, v, t) (11) 

t=0 t=0 

Furthermore, we notice that: 

Y Tp{f-v,v) - Y = ^^(S OU i,w E ) -^2jp(E in ,-VE) (12) 

reV,v feV,v S E 

Indeed, the first sum over the surface includes all the terms which are in the first sum over 
the volume, but not in the second, whereas it is the opposite for the second sum over the 
surface. This equation slightly simplifies: 

^lp(f-v,v)- Y Vi^v) = ^^(S OU i,w E ) - Sip(0) (13) 

f6V,t) rdV,v s 
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FIG. 2: The surface is totally reflecting; to compute the average return time to a site, we compute 
the average return time to the surface surrounding the site, and add 1. 

Finally, we can easily extend the model to the case where the surface portions may be 
either absorbing or reflecting. For example, it may be useful to compute the mean return 
time to a site, if the boundary conditions are totally reflexive, (see fig |2J) 

We denote S the points of the absorbing surface and £' the points of the reflexive surface, 
defined by the reflective boundary condition: 

pCfV.tlE^-tfc) =p(r* t^flE^tkO (14) 

The equation ^] is thus modified the following way: 

^2<p(f-v,v)- *Y^Tp{r,v) = ^^(S OU i,w E ) - ^^(E in , -fa) (15) 

r£V,v r£V,v E E 

E' E' 

The two last terms are exactly equal, and thus we have: 

^2<p(r-v,v)- Vi^v) = ^^(S OU i,{J E ) - S(p(0) (16) 

r£V,v r£V,v E 

where S is the surface of the absorbing portion. 

Thus, in the sequel of the article, we have to keep in mind that all the results apply to a 
lattice where a part of the boundary is reflecting. 



III. FIRST EXIT TIME 



We can now proceed to the computation of the moments of the first exit time. We have 
the following equation for the conditional probabilities: 

pif , v , t + l\r — v, v) = pif , v, t\f, v) H \p(f , v 1 , t\r, v") — pif 1 ', if, t\r, v) (17) 

v" 

This equation is simply the translation in terms of conditional probabilities of the rules for 
the behaviour of the random walker. We may at once sum over all the f and if on the 
boundary, as indicated in equation (j3J), and we have: 



pif- v, v, t + 1) = pif, v, t) H \p(f, v", t) - pif, v, t) 

v" 

We may now use the Laplace transforms: 



e s p(f— v, v, s) = p{f, v, s) + — Y^ \p(f, v" , s) — p{f, v, s 



When we sum over all f and v, the two last terms cancel out: 

^Tp(f,v") -^p(r,v) = 

and we can make use of the equation ()12j) to obtain: 



^]3(S ou4 ,{7 s , s) - y )p(Ein, -vs,s) 



[1 - e s ) ^2p(r,v,s) 



Thus, we have: 



S ((e- s % ~ 1) = (e- s ~ 1) °dV (e~% 
If we develop each side in powers of s, we obtain: 



fc=i 



U — ' mi 

1=1 m=0 



If we identify, we have: 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 



ill! MA = ZE^Y" I (t k - l \ 

k\ ^ ;s s j^i\(k-iy. \ I 



V 



(24) 
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And finally we have the relation: 

= a -f £ (,;) <'"- m >v (25) 

It is possible to obtain this expression directly from the evolution equation of t n . However, 
the computation is slightly longer, and we will not detail it here. 

An important consequence of this relation is that, in the special case of the first moment, 
we have the following simple result: 



<*> E = (26) 



S 

This explicit result is quite similar to the result obtained in the continuous case [16|: it 
has the same dependance on the surface and volume, but the numerical prefactor is modified. 
The simplicity of this equation makes it very easy to use. For instance, in the case of fig. 
El we obtain the result that the average return time is simply V. It depends neither on the 
frequency A nor on the shape of the volume! 

However, we must note that, for the higher-order moments, the results obtained are 
different from what we have in the continuous case : in the latter case, the nth moment 
of the surfacic first exit time depends only on the (n — l)-th moment of the volumic first 
exit time, whereas, on a lattice, we must take into account all the (n — 1) first moments of 
the volumic first exit time. 

Furthermore, we can obtain a lower bound for (t) v , if we inject the value of (t) s into the 
equation for the second moment: 

<a = ^(«v + i) (27) 

Since we have: 

<* 2 ) E >«^) 2 (28) 

we finally obtain the following bound: 

Wv >-^~\ (20) 

Note that there may be equality only in the case when every trajectory from surface to 
surface has the same length. For a square lattice, this is only the case for a ballistic motion 
(i.e. A = 0) when the volume has a square shape. 
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IV. RESIDENCE TIME IN A SUBVOLUME 

In many problems, the interesting quantity is not the time spent in the whole volume, 
but rather the time spent in a subpart of this volume. For instance, the random walker 
may be a chemical reactant which may react exclusively on several catalytic sites. Then, 
the subvolume corresponds to the ensemble of these catalytic sites. Here we will only derive 
the first moment of the residence time in this subvolume V, since the other moments are 
of limited interest, and have a much more complicated expression. What we will compute 
is the average residence time in V, assuming the random walk starts and finishes at the 
boundary of V. Thus, here, we will slightly modify the definitions of p(r,v,t) and t(f,v). 
The time will be the time spent in the subvolume V. We thus have: 




(30) 



if r G V . 




(31) 



if r g V 

Thus, we have: 




(32) 



if r e V . 




L" 



(33) 



if r £V 

We sum all this, which leads to: 



J^t(f- v,v) = ^2t(r,v) H ^2 [ t( f^) 




t(r, v)] + a D V 



(34) 



The equations (|16|) and (|20j) still apply here, and we have: 



(35) 



S 
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Thus, the residence time in a subvolume is proportionnal to its volume, and doesn't 
depend of the shape of the subvolume. We can also say that each site is visited in average 
^ff- times. We check that, if the subvolume is in fact the whole volume, we find the same 
result as previously (see equation (jZEj) ) 

V. NUMBER OF HITS ON A REFLECTING SURFACE 

Another question which may be asked is the following: On average, how many times does 
a particle touch a portion E" of the reflecting surface before exiting? If we note t the number 
of times the particle touch the surface E" before exiting, we have the following equations: 
(we note in this section E' the rest of the reflecting surface, and E is still the absorbing 
surface) 

p(r - v, v, t) = p(r, v,t) + — V [p(r, if, t) - p(f, v, t)] (36) 

V 

and the boundary conditions: 

p{X in ,-0 E ,t) = 6(t,O) (37) 

p(Kn,-vn>,t) =p(L' out ,vv,t) (38) 

p(El,-^,Ul)=p(C^,i) (39) 
Using the Laplace transforms, we get the following equations: 

p(f — v, v, s) = p(r, v, s) + — \p{r, if, s) — p(r, v, s)} (40) 

w 

p(E in ,-vx,s) = 1 (41) 

p(S-n) S ) = PiKut, ve>, s) (42) 

e s p(S-'„, -vz»,s) = p{Kuu vx»,s) (43) 
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Here we can sum the relation (J4*U|) over all f and v, and use the relations ()12p and (|2U|). 
which yields: 

- p(Sm, -^s)] + ]T [p(SL t , t&) - p(^ n , -ifc)] (44) 

S E' 

E" 

We define: 



6 E» 

Using the boundary conditions, we get: 

5 «e- st ) s - 1) + S" (1 - c -) (e-%, = (46) 

If we develop this expression into powers of s, and identify, we finally get the following 
relation for the moments of t: 

= f E (f) ( 47 ) 

m=l ^ ' 

We may notice that, for the first moment, we get the simple result: 

(*>e = f (48) 
This relation is again quite simple, and may easily be directly useful. 

VI. CASE OF A NON-UNIFORM ENERGY LANDSCAPE 

The most simple models of random walks, which we have considered until now, assume 
that all the points of the lattice are equivalent. But, in some applications, we have to take 
into account the fact that the vertices of the lattice may be at different potentials. For 
instance, in the case of a particle (vacancy, adatom) moving on a crystalline lattice, the 
presence of a inhomogeneity (typically, another kind of atom, which would be adsorbed in 
the surface) may modify the effective potential for our random walker around it. In the 
model we will introduce in this section, the various points of the lattice may have different 
energies. To take this particularity into account, we add to our model a reflexion probability. 
When going from the site x to the site y, the random walker has a certain probability R x ^ y 
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to be reflected, and a probability T x ^ y to be transmitted. The probabilities satisfy a detailed 
balance relation, which is: 



T p-E* — T P ~ E y — A 
- 1 x^y^ - 1 y^x^ x,y 



(49) 



(We condider kT — 1, and scale the energies accordingly) 
Thus, the law is: 



f(t + l) = f(t)+v(t) 

v(t) prob. 1 — A 
random prob. A 



v(t + 1) 



prob. Tp( t )^r(t)+v(t) 



(50) 



v(t + 1) 



f{t + 1) = f[t) 

-v(t) prob. 1 - A > prob. R^ t )->r(t)+v(t) 
random prob. A j 
As for the boundary conditions, there may be reflexions on the entrance of the volume. 
We note A-£ out ^ in = ^4s- Of course, if the particle is immediately reflected, the total time it 
will have spent inside the volume will be 0. 

We also redefine the average residence time: it will be the residence time weighted by the 
Boltzmann factors of the entry sites: 

E S e-^-t(S out ,t7 s ) 



<*>E = 



It is this average residence time we will compute here. We have: 



(51) 



p(f - v,v,t + 1) = 7>_£ 



+ (1-T, 

if f — v G V. As for the boundary conditions, we have: 



p(S out ,{7 E ,t + 1) 
p(S OU i,w s ,0) 



^7~v jV 

'-'out ,Zj tn 



+ ^7 V 1 , t) - p(S in , U E , *)] 



(52) 



(53) 
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p(£ in ,-<u s ,t+ 1) = Rz in -+x out 



~ Eir K S »> ^> *) - P( s *»> *)] / (54) 



Then, from this model, it can be shown (see Appendix [X} that: 



Er e " 



-Ef, 



^out 



This relation may quite easily be understood intuitively: If the energies in the volume are 
lower than the energies in the surface, the random walker will have more difficulty exiting 
the volume, and, thus, will stay longer inside the volume: the sum Ef e_Ef? w * n increase. 
Inversly, if the energies in the volume are higher than the energies in the surface, then the 
particles will tend to be reflected immediately, and the average time spent in the volume 
will decrease drastically. Moreover, if the energy landscape is flat, the energies are identical 
everywhere in the volume and the surface, we go back to the simple result Ij'ifijl . 



VII. JOINT RESIDENCE TIME FOR TWO PARTICLES 

We will now consider not one, but several particles diffusing independantly, and we wish 
to compute joint residence times. This is interesting in the case where we have several 
particles which must be on the same site to react [20]. For example, it may be two chemical 
reactants, or a vacancy and an adatom. If the two particles are strongly interacting, they will 
react as soon as they are in contact. The interesting quantity here is the meeting probability 
P. We can't compute it directly, but, given the average interaction time, we may have an 
upper bound: 

P <<t> (56) 

(it is interesting, since the average interaction time will genraly be much lower than 1.) 
Otherwise, if the two particles are weakly interacting (i.e. have a small probability p of 
reacting each time they meet), then the reaction probability will be approximatively 

P = p<t> (57) 

Thus, it is always interesting to compute the joint residence time. So, we consider here two 
particles, which may either start from the bulk or the boundary of the volume. We define 
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the joint residence time as the amount of time spent by the two particles on the same site 
before one of them exits. Thus, we can define q(f,v,f',if,t), the probability that the two 
particles meet t times before they exit, given that the first one starts from the position r 
with the speed v, and the second one starts from the position f with the speed if. We can 
also define the average interaction time given the initial positions and speeds t(r,v,r , ,v f ). 
We have the following equation: 



(The 5(r, r*) is the Kronecker delta function, whose value 1 if r — f and otherwise) 
If we sum this over all f, r* , v, v 1 ', we notice that the right-hand side of this equation nicely 
simplifies: 



The last term is simply the number of quadruplets [r,r',v,v r ) where r = r* '. We have the 
following relations, similar to equation (J16|) : 



y^Jjr- v,v, 7* - if, if) = ^ t(r, v,r'-if,if) + } j t(E 0Ut , r 1 - if, if) (60) 

+ y^ y t(E in , -v^, r 1 - if, if) 



s 

The last term in this equation is zero, since if a particle is at the position Ej„ with the speed 
—vj], it immediately exits. We thus have: 




(58) 




(59) 




(61) 




Finally, reporting this result in the equation (JBUj) . we have: 



2 Y t(Z outl vx,r ,if) + ^£(£ oui ,{T s ,£' out ,^ E ) = a 2 D V 



(62) 
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If we define the following average joint occupation times: 

(*>sv = TF~~q Yl t(Z°u*,VEy,tf) (63) 



Va D S 

is the average joint occupation time when a particle starts from the volume and the other 
from the surface; 

is the average joint occupation time when the two particles start from the surface. Then, 
we have the following result: 

a 2 D V = 2a D VS (t) sv + S 2 (t) sa (65) 

Thus, we cannot have a totally explicit result for the joint residence time; we just derived 
a relation between two different joint residence times, depending on where the two particles 
start from. However, this relation gives immediately an useful upper bound to these joint 
residence times: 

(t) E2 < ^ (66) 

<*>ev < % (67) 

Since, in trial cases, we found out that the magnitude of the two terms in the relation (}65|) 
was similar, we have good hope that these relations will at least give an upper bound which 
is of the good order of magnitude. 

It is possible to have a similar result for n particles: it is a relation between n different 
averages (depending on the number of particles which start from the boundary, and the 
number of particles which start from the volume. The result (see Appendix|B]for details on 
the computation) is the following: 

a D V = E (") S m {a D V) n - m {t)^ myn . m (68) 



where (t) Sm y„_ m is the average joint residence time for n particles, of which m start from 
the surface and n — m from the boundary. 
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VIII. CONCLUSION 

The results obtained in this paper significantly extend the results previously derived 
Q 3- They show that the mean first exit time behaves differently for a discrete lattice 
and for a continuous media, not only quantitatively (as for equation (J2HJ), for instance), but 
also qualitatively (see the relation for higher-order moments (|25Jl ) It should be pointed out 
that we obtained explicit exact results (equations (|2B |) .(|55 |) .pS ]) and which is not so 

common in the bounded random walks literature. Furthermore, they apply to all kinds of 
lattices and random walks: because they have a very general range of application, these 
results can be very useful when the evolution equations cannot be solved exactly. In fact, it 
is even possible to generalize our methods to irregular networks, where the connectivity can 
vary from one site to another. This may be of special interest for complex networks which 
can be found in ethology , economy, neural networks or social sciences. Such an extension 
is in progress. 

APPENDIX A: COMPUTATION OF THE MODEL WITH ENERGIES 

We can rewrite the equations (|53|) and in terms of average time: 

t{f -v,v)-l = (l{f, v) + ^~Y^ P( f ' ^ _ ^ r "' ^] j ( Al ) 




if r — v G V, If we define: 

a s = £(£ in , vx) + — J2 [t(Sin, &) ~ *(£*«, ^e)] + 1 (A2) 

Then we have: 

t(Z out , vz) = T Sout _ Sin Q; s (A3) 



{(^-^ = (1-^0^ (A4) 
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If we sum the relation (jAl|) over all r and v, weighted by the appropriate Boltzmann factors, 
we get: 



r,v,r—v£V 



E A r,f-v \t (f, V) - t(f - V, -V) E 



r,v,r—vdV 

+ E e~ Ep -" [ t(f - v, —v) 

r,v,r—v£V 

We will see how all these terms may transform. 



t(f, v') — t(f — v, v') 
+t(r — v, —v) — t(f, v) 

+ — \t(f — v, £>*) — t(f — v, —vj] ] 

E v' / 



E e- Ef -n(f- v, v) = e~ E %r, v) - E e _i?E ™t(S m , -v s ) 

r,v,f—v£V f,v £ 

We also have: 



f ,v,f— v(zV 



E A??-jt{f- V, -V) = E v') - E ti^in, 



r,v,r—v£V 



(AS) 



(A6) 



(A7) 



(A8) 



(We just take v' = —v, f = r—v). Note that these two expressions are identical only because 
of the detailed balance condition, which implies that the function A has to be symmetrical. 
It is important since these two terms must cancel each other. We must also compute: 



} J A^_ ff t(r, if) = } u A^fft(r, jf)-/ u 



(A9) 



r,v,r— «gV 



E Ap?-ijb(r- v,v>) = 22 Ap"-ff> ti i"t(r",v') - E*( s fa>^) (A10) 

(We just take v" = —v, f" = r — v) These two terms also exactly cancel each other because of 
the detailed balance condition. We have of course similar relations with the terms in e~ Eit - s 
We now use all these relations in our main equation, which yields: 

E e~ K w, v) - E e ~ ESi <z in , -tfe) - e e ~ Er + E e ~ ESm = 



r,v 



A 



V e~ E H{r, v)-Y] e- E ^t(^ in , «e) + — E e ~ E *W, ~ W, v)\ 

r,v 1j r,v,v' 

+ ^EE e ~ ESm P( S *>» ^) - *)] 



(AH) 
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We use the relation (jA4j) on the left-hand side, and recognize «s on the right-hand side, we 
get: 

E Ffi E E 

Since, because of the relation (|A3J) . we have 



We - ■ 

We get the relation: 



(*>£ = ~ X ( Al3 ) 



(*>E = ^fi^ ( A14 ) 
Z^E e 

APPENDIX B: COMPUTATION OF THE JOINT RESIDENCE TIMES FOR N 
PARTICLES 

Here we will evaluate the joint residence time for n particles, i.e the amount of time 
they will spend all in the same site. We will not write the equations in the fullest detail, 
since it would take pages, but the calculation is quite similar to the calculation for only 
two particles. We will define the times t(fx, V\, . . . , f n , v n ) which are the average time the n 
particles will spend together, given that the kth particle starts at the position 7% with the 
speed Vk- These times obey a relation quite similar to the equation (J58|) . but this relation 
is quite unwritable in the case of n particles. However, we can see that the right-hand term 
will simplify the same way it does for two particles, and we will get the equivalent of the 
relation (fo^j) 

Y Kn-v 1 ,v 1 ,...,r ri -v n ,v n )= Y t{n,v 1 ,...,f n ,v n )+a^ ) V (Bl) 

fi,vi,...,f n ,v n fi,vi,...,r n ,v n 

Again, a^V is simply the number of 2n-uplets (fx, vx, ■ ■ ■ ,f n , v n ) which satisfy fx = ■ ■ ■ — f n 
We will also have the relation, similar to the relation (|61|). where we put together the terms 
which are similar: 

Y t(fx-vx,vx,...,f n -v n ,v n ) = Y t(n,Vx,...,r n ,v n ) (B2) 

fl,vi,...,f n ,v n fi,vi,...,r„,v„ 

sr^ f n \ tyvW -CO v( m ) -^ m ) - - - \ 

"t~ / J I I / J / l \ 2 ^out) V T, J • • • ) ^out ) V E > r m+l, Vm+X, ■ ■ • > r n, Vn) 

m=X ^ ' EC 1 ),...,E( m ) r m +i,Vm,+l,-,r n ,Vn 
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The average joint residence times are defined by: 




(B3) 



We thus get the relation: 



a D V = J2( )S m (<J D V) n - m {t) 




(B4) 
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